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Chapter 18. Gordon decomposition of the vector/axial currents 



18.1 Phase gradient versus Magnitude gradient 

Wave-functions are constructed from complex amplitudes with both a 
phase 4> and an absolute amplitude a. Both components play different 
roles which we want to study in more detail. We want to separate certain 
quantities into phase related and magnitude related components. Impor- 
tant in this respect will be the Gordon decomposition which separates the 
vector current in a charge and a spin related component. 

In this section we use a more elementary approach using the Klein Gordon 
equation. Each component of the Dirac/Weyl spinor obeys separately also 
the Klein Gordon equation. To show this first we start of with the Dirac 
equation. 

A A ( iPl \ _ ( ipL 



We now square the operator on the left hand side. In the first order 
equation the tpL and -ipR chiral components are each the source of the 
other coupled via the constant m. Applying the operator twice couples tpL 
and V'i? back into an expression where they each are their own source. 



)d 
cr" l U » 



The result is that ipL and if)R are given by independent equations. From 
the commutation rules of the Pauli matrices we derive. 

a* \ ( a- \ f / if /i = i/ 

or" M er" " 1 if fi ^ v 



Which shows that all four components obey the Klein Gordon equation. 

Dip 1= -m 2 ij I (18.4) 



This makes if physically meaningful to study the relating between the phase 
gradient (j)^ and the gradient of the absolute amplitude a in the context of 
the simpler Klein Gordon equation first. 
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We can describe wave-function around an arbitrary point which we choose 
as the origin of a coordinate system x^ as an exponential with phase change 
rate 4>^ and absolute amplitude change rate a M . where 4>^ and a M are four- 
vector functions depending on x M . 

$ = eK+^K (i8.5) 



The Klein Gordon equation becomes. 

( 8^ + m 2 ) e K+^K = o (18.6) 



We evaluate this and then set = because we are interested in this 
particular point only. This gives us the following expression. 

(a^ - (f>^ + 2i - 2i &* (a M + i<f>„) ^ if) = (18.7) 



We can spit this expression into two independent equations, one for the 
real and one for the imaginary part. 

-m 2 = + 2 real part 

d^cj)^ = —a^(j)^ imaginary part 

If a M = then this simplifies to the following two equations. 

— m 2 = Momentum space Klein Gordon eq. , , 

= Continuity relation 



The first one is simply the Klein Gordon equation in momentum space. It 
is however also locally true, that is, energy and momentum are allowed to 
be locally different as long as the equation holds in each point. The second 
equation is the continuity equation for energy /momentum or equivalently, 
the charge/current density which transforms in the same way. 
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18.2 Gordon decomposition of the vector current 

We now precede with the Gordon decomposition of the vector current of the 
Dirac equation. It splits the vector current Jy into a current related to the 
charge and a current related to the spin. The vector current, normalized 
as a charge/ density current is. 

jy = ~- (18.10) 

m 



Where the sign of e determines the sign of the electric charge. We have 
divided by mc and then multiplied by ec to go from the momentum four- 
vector to the charge/current density. 

We will split this current into one ( j£) which depends on the phase change 
rates but which is independent of the changes of the magnitude, and an- 
other (J a) which does not depend on the phase change rates but only on 
the changes (derivatives) of the magnitude. 



This is done by either adding or subtracting a complex conjugate term. 
<ft = ^( (a + i<f>) - (a + icft)* ) 



a -- 2 



I (a + i(f>) + {a + i<P)* ^ 



(18.12) 



We achieve this by first splitting the vector current into two equal parts. 

j v = -JL^yv + ) (18.13) 



We replace ip in the first term above with the help of the Dirac equa- 
tion, and we replace (ft with the use of the complex conjugate of the Dirac 
equation. The Dirac equation and its complex conjugate are. 



i^dyift = mcip, 



— irfdvift = mc ip 



(18.14) 
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The 7^ matrix is just a constant for differentiation so we can move it 
through the d v operator, we get. 

■ji ieh f 7 dip d$ v \ 

J y = ~2^c { ^ 7 ~ a^ 7 7 ^ J (18 ' 15) 



These two terms now subtract if fi = v and they add if jjl / v because. 

Y l u = l v ' 7 M if At = ^ MSIfil 

7 m 7 ^ = -Y 7^ if /x / i/ 1 J 



The current related to the phase change rate becomes (using 7°7° = 1 and 

7 Y = -l) 



We recognize this as the expression for the charge/current density of the 
Klein Gordon equation. Subsequently, the current which depends on the 
magnitude changes is. 

« = -^^(^) (18 - 18 » 



Where we have used the definition a^ v = | (7^7^ — 7^7^)- The deriva- 
tive here is acting on the whole bilinear t^a^ip, from which the phase 
information is already eliminated. 

We showed in section ?? that the bilinear tjja^ u ip transforms as the elec- 
tromagnetic field tensor, or equivalently, the polarization/magnetization 
tensor. Indeed, when we recall Maxwell's inhomogeneous equations in ten- 
sor form, (with e Q , fi and c set to 1) 



3<t> 



-d v F^ 








—E x 


-Ey 


-E z 
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E x 





-B z 
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dx v 


Ey 


B z 





-B x 




Ez 


-By 


B x 






(18.19) 
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We see that the divergence of the rows does lead to a charge/current den- 
sity. This means that we can define tpa^^ip as an electromagnetic field 
tensor tensor. In the rest frame ipa utJ "ip becomes. 



f3=o 



o 
o 









-is, 





is z 




tSqi % S 7 





-ISy 

IS X 





(18.20) 



Where s is the unit spin vector. Thus, there is only a magnetic field in the 
rest-frame. Under an arbitrary boost [3 the tensor becomes. 



~2~ 





i Pi )x 

i Pi )y 

>Pl )z 



- ( /?7 ) x 


-i(s\\ + s±~/) z 
i(h +s±l) y 



" ( s® Pi )y 
i{s\\ + s_l7) 2 


-i{s\\ + s±<y) x 



- ( s /?7 ) z 
-i{s\\ +s±j) y 
i(s\\ + s±7)x 


(18.21 



The first column and row now contain the E components which correspond 
to the transformed B field components. We note that the B field compo- 
nents are imaginary while the E field components are real. This leads us 
to define the electromagnetic field using complex numbers: 

F = E + iB, or F = E + B (18.22) 



That is, we have to multiple the B fields with an extra factor i in order to 
obtain the correct electromagnetic tensor. 
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(18.23) 



This expression corresponds with the usual inhomogeneous Maxwell equa- 
tions given in SI by. 



p = e Q V • E, 
p = V-D, 
P = VP, 



Ja 
Ja 
ja 



1 „ ft #E 

— V X B £o "777 

Mo 



V x H - 

V x M - 



(18.24) 
(18.25) 
(18.26) 



Using either the E and B fields or the D and H fields. We will however use 
the polarization vector P and the magnetization vector M as the physical 
interpretation, where the spin is considered to give rise to the magneti- 
zation of the wave function. M is strictly zero outside the wave-function 
unlike H and B. 



18.3 Effective spin-current around the wave-packet 

Figure 18.1 shows the axial current of a Gaussian like shaped wave packet 
at the left and the corresponding effective current at the right. Recalling 
Stokes law we see the correspondence between the two. 

In the center, the axial circular currents cancel each other out. At the 
edge, where there is a gradient, an effective current arises as a result. This 
effect also occurs in magnetic media where an effective current occurs due 
to the change in magnetization, j = V x M 

The quantity -ipa utl if) is therefor best interpreted as the (relativistically 
transformed) magnetization of the electron's wave-function. Note that 
tpa v,J "tp is zero if there is no local charge density since it is proportional to 
the absolute square of the amplitude of the wave function. 
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These insights are due to Belinfante [?] (1939), JJ.Sakurai [?] (1967) and 
H.Ohanian [?] (1986). 




Figure 18.1: Axial and effective current density of a Gaussian packet 



From equation (18.20) we see that, in the rest-frame, the axial current 
= ?/>7^V is proportional to the magnetization components of the tensor 
field vpa UfM ip. The two however transform differently under a boost. 

Total magnetic moment: ft = Jm dS = \ J rXj dS (18.27) 

We can interpret the left side of figure 18.1 as the magnetization field M 
and the right side as the effective current j=V x M. (In the rest-frame). 
The total magnetic moment jl can be derived from both as shown. 



18.4 Charge/magnetic-moment ratio of leptons 
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18.4 Charge/magnetic-moment ratio of leptons 

The table shows the mass/spin ratios and the charge/magnetic-moment 
ratios for the three leptons. The electron, the muon and the tau-lepton. 
The ratios for the electron are normalized as 1:1. 

These ratios are explained by the Gordon decomposition of the Dirac vector 
field, except for the small magnetic moment anomaly terms caused by 
interaction with virtual particle pairs. 



The ratios correspond to the difference between the two currents j£ and 
ja- The first current does depend on the phase change rates while the 
second current doesn't. 



leptons: 


ELECTRON 


MUON 


TAU 


Spin z-component 
Mass in MeV 
Mass ratio 


1/2 h 

0.510998918 
1.000000000 


1/2 h 

105.6583692 
206.7682838 


1/2 h 

1776.99 

3477.48 


Charge 

Magnetic moment 
Mag. moment ratio 
Magnetic anomaly 


1.0 

0.928476412c-23 

1.00000000000000 

1.00115965218085 


1.0 

4.49044799e-26 
1/206.7669894 
1.00116592080 


1.0 

2.6700053e-27 

1/3477.43 

1.00117324 1 



j£ corresponds with the energy /momentum vector as well as the charge/ 
current density vector, both transform in the same way. 



Being dependent on the phase change rate of the wave functions means 
that they are larger by a factor proportional to the mass in the rest frame 
compared to the quantities which correspond to The spin and the 
magnetic moment density. 



1 theoretical value from QED and QFT 
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18.5 Gordon decomposition of the axial current 

Similar to the Gordon decomposition of the vector current, we will now 
handle the Gordon decomposition of the axial, current. The axial current 
is given by. 

3A = ~ VVW (18.28) 



Where the sign of e determines the sign of the charge. We will again split 
this current into one (jA<j>) which depends on the phase change rates but 
which is independent of the changes of the magnitude, and another {jAa) 
which does not depend on the phase change rates but only on the changes 
(derivatives) of the magnitude. 

3A = 3A4> + JAa (18.29) 

The first step is to split the axial current into two equal parts, 

3 A = + ^VtV ) (18-30) 



and then to replace replace ip in the first term above with the help of 
the Dirac equation and to replace <j) with the use of the complex conju- 
gate of the Dirac equation. The gamma matrices are just a constants for 
differentiation so we can move them through the d u operator, we get. 

* - -*k ( M t ~ w^-r* ) < 18 ' 31 ' 

These two terms now add if (i = v and they subtract if fx ^ v because. 

7 M 7V = -T^V if/i = i/ , , 



The signs are exactly the opposite of what they are in the case of the 
Gordon decomposition of the vector current, due to the way 7 5 operates. 
Now, the current which does only depend on the phase change rates and 
not on the magnitude change rates becomes. 



18.5 Gordon decomposition of the axial current 
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JA<p 



(18.33) 



Where a^ 5u is the dual tensor of a^ u and p u represents the phase change 
rates. Note that p v is covariant now Pu = ( P t , ~P X , ~P y , ~P Z ) while it is 
contravariant in case of the vector current where the spatial components 
obtained an extra minus sign because of 7*7* = — 1 while 7°7° = 1 



The part of the axial current which depends only on the changes of the 
magnitude and not on the phase change rates is. 



eh 







jAa 



2m 2 c dx v 



■07 ip 



(18.34) 



Where f/>7 5, is zero under quite general conditions for a free electron. 
This leaves as jA0 as the basic contributor to the axial current j'a- We 
have defined the dual tensor a^ 5u as. 



a ^u = * [ yy^i/] ? with ^5,Y] = {^7 5 Y -7 U ^7 5 ) (18-35) 



or, written out in full for your convenience: 
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-ia z 



1(7 
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-a z 
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-a* 
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ia z 



-rrV 









<7 X 



(18.36) 
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The bilinear field V> a^ v tp relates to the bilinear field ip a^ v tp in the same 
way as the electromagnetic field tensor F relates to its dual, that is, the 
relationship, 

^a^ip <=> ixr^ij (18.37) 



is equivalent to the relationship between F^ v and *F^ U : 
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(18.38) 
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- D 2 
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T» = ^— F» a A a (18.41) 

47T 2 

3a<$> = 3u + (18.42) 

We recognize this as the expression for the charge/current density of the 
Klein Gordon equation. Subsequently, the current which depends on the 
magnitude changes is. 
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